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SUMMARY 

An expression is developed for the lift of a slender afterbody in 
terms of the "slender-body" approximate potential at the after end in 
conjunction with a suitably calculated value of the potential at the 
forward end. The failure of the usual "slender-body" theory to predict 
any lift on a slender cylindrical afterbody is thereby corrected. The 
same expression is used to compute the part of the interference lift 
generated by the presence of a neighboring body, due to the interference 
upwash. Another expression is developed to compute the remainder of 
the interference lift; due to an interference pressure gradient. The 
lift is determined for a cone- cylinder body; a cylindrical afterbody 
of a slender wing-body combination; and three combinations of bodies 
to illustrate the method. 


INTRODUCTION 

In the search for better aerodynamic configurations for supersonic 
flight; the missile shape has become more complex. One of the problems 
arising is the effect of interference between bodies which occurs; fox 
example; when the propulsive unit is mounted external to the fuselage. 

An adaptation of slender-body theory for obtaining a simple approxima- 
tion of the interference lift between such bodies was developed at 
the NACA Lewis laboratory and is presented in this report. This approxi- 
mate theory is also applied for the calculation of the lift of an afterbody, 
which is herein defined as the portion of the body behind the wing trailing 
edge for a wing-body combination or as the portion of the body aft of the 
nose section for a wingless body. 

Slender-body theory as originally developed by Munk in studying 
the lift of airships (reference l) has proved useful in predicting the 
lift of low-aspect-ratio wings, slender bodies, and slender wing-body 
combinations at supersonic speeds (references 2, 3, and 4). However, 
slender-body theory yields the unrealistic result that a cylindrical 
afterbody of a wing-body combination and a cylindrical afterbody of 
a body carry no lift. 

Although slender-body theory may be grossly in error on the lift 
of a slender afterbody, the theory closely approximates the correct 
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value of the part of the surface potential at the rear of the afterbody 
proportional to the angle of attack. An expression utilizing this fact 
is developed herein for the afterbody lift in terms of the slender-body 
potential at the after end in conjunction with a suitably calculated 
value of the potential at the forward end. The more accurate methods of 
calculation of greater labor (for example, linearized theory) are thereby 
limited to the forward or nonslender portion of the body. 

The same expression is used to compute the part of the interference 
lift of bodies of revolution due to the upwash field generated by each 
body on the other. Another interference effect considered herein is 
that due to the variable pressure field generated by one body in the 
vicinity of the other, which results in a buoyant force. 

As examples of the method, the lift of a cone-cylinder body, the 
lift of a cylindrical afterbody of a slender wing-body combination, 
and the lift of three combinations of bodies are determined. 


ANALYSIS 
Basic Equations 

Consider a slender body of revolution at zero angle of attack in 
a uniform supersonic stream of velocity U. There will be a certain 
disturbance velocity field and an associated velocity potential produced 
by the body. Because of symmetry, this velocity field will give rise 
to no lift. 

Now let the body be given a small angle of attack a and be subject, 
in addition, to a small disturbance field (u,v,w) generated by a second 
body (fig. l) . (The symbols used herein are defined in appendix A. ) 

The first body will develop an additional velocity potential due to a 
and to the disturbance field. If v and w are suitably restricted 
it is shown in appendix B that to a certain accuracy (a) only that part 
of this additional potential due to a and to w will give rise to 
lift, and (b) cross product terms involving both potentials may be 
neglected for the lift. Call this lift-producing incremental potential <P. 

Equations of motion and boundary condition . - The incremental 
potential <p is presumed to satisfy the governing equation for small - 
disturbance, isentropic irrotational flow 


P ^xx “ ^yy “ ^zz " ® (■*■) 

In addition, <p is specified to satisfy the boundary condition that the 
cross flow, Ua +. w, must be canceled at the surface of the body; that is, 
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(‘Pr)^ = - (Ua + w)sin 0 ' (2) 

on the body. This boundary condition is linearized, not exact. 

Pressure coefficient . - That part 6p of the entire pressure 
increment Ap which contributes to the lift is, to the present accuracy, 
given by the linearized Bernoulli equation 

8p _ -Z(V X + u) 
q IT 

Specifically, the accuracy is such that all terms of order a dR/dx or 
(dR/dx) but no higher, which can contribute to the lift will be 
included; this result is proved in appendix B. 

Lift . - The axis of the slender body under consideration will 
always be taken parallel to the x-axis. For simplicity the sidewash v 
and w generated by the interfering body are considered to vary only 
with the coordinate x in the general vicinity of the primary body; 
this restriction is the one mentioned earlier. (The sidewash will not 
be considered further because it contributes nothing to the lift 
(see appendix B) ) . The approximation should be valid if the bodies 
are not too close to each other. 


The lift of the body may be obtained by integrating the component 
of the pressure in the lift direction around and along the body. The 
expression for the lift is 



where R is the radius of the body and x^_ and xg are the x- coordinates 
at the start and base of the body, respectively. 

Substituting equation (3) into equation (4) gives 



The second integral in equation (5) represents the lift of a body 
immersed in a variable pressure field and may be considered as a buoyant 
force. The integral will be discussed in detail in the section on 
interference between bodies in combination. 


( 3 ) 
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Development based on first Integral in lift equation . - This 
section is concerned -with the transformation of the first integral in 
equation (5) into a more useful form involving the surface potential. 
Hereinafter this first integral -will he called L. 

For constant 0, the quantity Ru dx may he written 


Ru dx a R dj> - R(<p dr) = d(Rcp) - <p — dx - R(<P r ) ~ dx 


dR 


r»R 


dx 


r=R dx 


( 6 ) 


Insertion of eqmtion (6) into the first integral of equation ( 5 ) yields 

f f ^ 9 { B ^ - - /*[• - Jf «*] 


d0 (7) 


0 


The equation for the upwash perturbation velocity is 


v =* (p r sin 0 + — <Pq cos 


( 8 ) 


Substitution of the value of <p r sin 0 from this equation into equa- 
tion (7) results, upon simplification, in the following equation for L: 


2jc x 2 2jt 

L 3 ~ J (*2*2 “ h 1^l) sin e ae - & J 1 dx Wd0 + 

0 x ± 0 


x 2 2jt 

25f Sfcf iC i . -a s 9 1 SB 

U J dx U d0 


0 


which reduces to 
2 * 


2it x 2 2jt 

L » ~ f* (^2 9 2 ~ R l v l) s in 0 d0 - ~ R ^ dx J* w d0 (9) 
0 x-i 0 


because the integral in the last term is zero. The first integral in 
equation (9) may he integrated directly because, as shown in appendix B, 
the potential is of the form 


<p = sin 0 G(x,r) 



MCA TN 2669 


5 


Equation (9) then becomes 

Xo 2it 

/ "I **X CIO) 


This equation is fundamental for subsequent developments . 

If <J>2. and <Pg> Ike values of surface potential at the forward and 
after ends,, respectively, of the section of the body under consideration, 
are known as well as the vertical component w of the body perturbation 
velocity, the lift may be obtained from equation (lO) . In the following 
sections, procedures are developed for determining the lift for cases 
in which the integration of the second term in equation (10) is negligi- 
ble or zero. ('When the second term is not zero, it is probably simpler 
to obtain the lift from equation (5) because the same type of computation 
is required to determine the w- distribution as is needed to determine the 
u-distribution. ) For these cases, <P2 is assumed to be the potential given 
by slender-body theory. ~ 

Slender-body theory for variable cross flow . - The central approxi- 
mation of slender-body theory is the following: In any plane x « con- 

stant, the value of the potential increment <P associated with the 
cross flow is the two-dimensional potential of a cylinder having the 
local radius and moving with the local cross-stream velocity. For 
supersonic flow with uniform cross flow (that is, a fixed angle of 
attack and no variable interference upwash), this result is derived in 
reference 3 by taking the limit of the linearized value for radius 
approaching zero. The corresponding result will now be derived in 
similar fashion for the more general case of variable cross flow. 


The potential in integral form mny be written 


x-f3r 

sin 9 P g(g)(x - g)dg 

0 ^(x - £) 2 - 0 2 r2 


(H) 


The expression may be interpreted as the potential of a distribution 


of doublets of strength 
appendix B). 


J 1 g(x)dx per unit length (g(x) = gp*(x) 


of 


When the distribution function for the limiting case R 0 is 
related to the local cross flow as in reference 3, there results 



( 12 ) 


T 
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Inserting tMs value of g(x) into equation (ll) for the potential, 
applying the slender-body approximation (radius approaches zero), and 
integrating yield 

<p => |jja + w(x)] - ■ - (13) 

On the hody (r = R), equation (13) reduces to 


<p = Jua + w(xf| R sin 0 


(14) 


Equations (13) and (14) represent the slender-hody approximation for 
the case of a variable cross flow w(x) . 

_2it 


Demonstration that 


f 

0 


v d0 = 0 for slender hody section. - 


The lift equation (lo) is particularly simple when the second term 

reduces to zero. This reduction occurs when either R » constant (the 

2« v 

section under consideration is cylindrical) or J"* w d0 » 0. It will 

2m 0 

now he shown that J"* w d0 = 0 is obtained when the afterbody is suf- 

0 

ficiently slender that slender-hody theory may he applied. 

The upwash may he obtained by inserting equation (13) into equa- 
tion (8) and performing the indicated differentiation. The result is 


w ta [jja + w(xjj R 2 /r 2 (cos 2 0 - sin 2 0) 
or on the hody (r » R), 


(15) 


w — [jja + w(xjj (cos 2 0 - sin 2 0) 


which immediately leads to the result 

2 it 


f 


w d0 ss 0 


( 16 ) 


wherever slender-hody theory is applicable. 
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Applications 

Slender pointed tody of revolution . - If, in addition to "being 
slender everywhere, the "body of revolution tapers to a point at the 

2it 

forward end, then w d0 = 0 and cp^ = 0. Thus the lift eq.ua- 

0 

tion (10) reduces to 

_ 2ixqR2<P2 

~ U sin 0 

With the slender-hody valtie for cpg (equation (14) ), the final result is 



L = 2qAg 



(18) 


where Ag is the area of the "base of the "body. Thus the lift of a 
pointed slender "body in an up, -rash field depends on only the upwash and 
the cross-sectional area at the base. 


Open-nose bodies and cylindrical afterbodies of bodies . - For a 
circular cylinder (cLR/dx = 0), the lift from equation (10) becomes 


2qatR 92 ~ <Pl 
U sin 9 


(19) 


If the potential at the base of the cylinder is given by the slender- 
body value (equation (14=) with w =■ 0), the lift is the same as that of 
a slender-pointed body with the same base area because is zero. 

Indication of the extent to which use of the slender-body value is 
justified may be obtained from reference 5 where an exact solution of 
the linearized differential equation ( equation (l) ) has been obtained 
using operational methods for the flow past finite bodies. Ward's 
results for the lift of a cylinder of length Iq indicate that the 

slender-body lift is realized within 2 percent for l 0 /W> 8. 

The lift of a cylindrical afterbody of a body is given by equa- 
tion (19) where <P^ is the potential on the body at the start of the 
aftersection. The lift of the afterbody would be zero if <P^ and 
<pg were both tahen to be the slender-body values because <Pj. "would 
then equal <P g- To obtain a more accurate value for the lift of the 
afterbody, tp-^ should be obtained by more exact methods, such as the 
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stepwise procedure of reference 3. In the EXAMPLES section, the lift 
of a cone- c ylinde r model is computed as a function of the cone angle hy 
this procedure. 

If in the vicinity of the lip the body is assumed to be quasi- 
cylindrical (dR/dx * 0), the second term in equation (10) will be 
approximately zeroj if the body is slender, the term will still be zero 

2 it 

farther downstream because I w d0 = 0 as shown. Lighthill points 

4 ) 

out in reference 6 that the discontinuity at the lip of an open-nose 
body affects only the region directly behind it for a distance of the 
order of the diameter of the body. The flow then behaves as if the 
open- nose body were a pointed body with a continuous slope. By use 
of operational methods, reference 6 shows that the lift of an open-nose 
body of sufficiently large fineness ratio is given by the slender-body 
value. 

Afterbodies of wing-body combinations . - The analysis used for a 
body of revolution may be considered to apply equally well to a wing- 
body combination provided the potential <p is no longer restricted 
to the form sin 6 G(x,r) . Consequently all the earlier formulas not 
depending on this restriction still apply. In particular, equation (9) 
must now be used in place of equation (10) . The lift of the portion 
of the body behind the wing trailing edge of a wing-body combination 
may be determined from equation (9) with <p^ equal to the potential 

on the body at the wing trailing edge. In order to satisfy the basic 
assumption of constant interference downwash and sidewash in each 
transverse plane around the body (due in this case to the wing and its 
wake), the diameter of the afterbody should be small compared with the 
wing span. Reference 7 gives the complex stream potential for a slender 
body of revolution with small-aspect-ratio triangular wings mounted on 
the cylindrical part of the body. The velocity potential is obtained 
by evaluating the complex stream potential on the body since the stream 
function is zero on the body. The result is 



where b is the maxi mum semispan of the wing and R is the radius of 
the cylindrical part of the body. The slender-body approximation to 
the surface potential at the base of the afterbody is (equation (14) ) 
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where v 2 is the induced upwash velocity generated hy the wing. These 
values of and cpg may now he substituted into the first term of 

equation (9) for the lift. The second term in equation (9) is zero 


•because 



w d0 vanishes as a consequence of the slenderness. The 


0 


result for the lift of an afterbody of a slender wing-body combination 
is found to be 


L = 2qa 

A 1 + -^2 ( 

1 + ^2 \ _ 2 R 1 (b 


Rl 2> 

o + -i- 

2 

) tan - - 1 - ( — J 


j. u y 

v UaV \ 

b / \ 

b / 

' \b/_ 


( 21 ) 


The lift of a cylindrical afterbody for the limiting case where 
the vortex sheet is displaced sufficiently far from the afterbody 
that the downwash induced by the wing at the base of the afterbody is 
negligible (w 2 /0a, * 0) is compared in figure 2 with the lift of a 
slender wing-body combination with no afterbody. The assumption of 
small R/b is violated in the region where the radius of the body 
approaches the semispan. However, the lift of the afterbody approaches 
zero in this region and since the lift given by equation (21) also goes 
to zero, equation (21) may be considered a good approximation to the 
lift over the complete range. The lift of the wing-body combination 
is obtained from reference 4 and is presented as the ratio of the lift 
of the combination to the lift of a pointed low-aspect-ratio wing. The 
lift of the wing is 

%, - ^ ( 22 ) 

Combinations of bodies . - Because the linearized differential 
equation of motion is assumed to apply, the solution for a combination 
of bodies is a linear superposition of the solutions for the bodies 
alone and for the interference effects between the bodies. The boundary 
condition to be satisfied is that the normal velocity at the surface 
of the bodies is zero. 

In order to determine the interference lift for the bodies in 
combination it is necessary to know the disturbance one body produces 
on the other. In particular, the upwash and pressure field of one 
body contribute to the lift of the other body. The interference pres- 
sure field can be computed from the interference upwash field by 
Integration of the irrotationality condition 

chi dw 
Sz cbc 

to determine uj the pressure is obtained from equation (3). 
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Computation of interference upwash 


Several methods may "be used to calculate tlie upwash field generated 
"by a tody at zero angle of attack and at angles of attack. For a 
tody at zero angle of attack the method of reference 8 may te applied 
to determine the upwash. The velocity potential for the supersonic 
flow atout a tody at zero angle of attack is given ty 


X 



f(Qag 

^ (x - £) 2 -8 2 r 2 


(23) 


where f (x) is the local strength of the source distribution. The 
integral expression for the upwash is 


w 


x-pr 

sin 8 

0 


f'(g)(x - g)dg_ 
^(x - g) 2 - p 2 r 2 


(24) 


A first approximation for the source strength of a slender tody is 


f (x) = UR (25) 

dx 

Equation (24), with the source strength given ty equation (25), can te 
integrated directly for a given tody. 

An analogous procedure may also te applied for nonzero angle of 
attack; that is, the slender-tody value of g(x) (equation (12) with 
w s* 0) for a given tody is substituted into the integral expression for 
the upwash 

w _ sin 2 9 p ^ g*(g)(x - g) 2 dg + cos 2 9 p ^ g(€)(x - €)dg 

r2 J o 'Jc* - « 2 - gV J o Vc- - «) 2 - 

(26) 

and the integration is performed without letting the radius approach 
zero as before . The result is thus dependent on Mach number, whereas 
the upwash computed ty the slender-tody approximation (equation (15) 
with w = 0) is not. 

A more accurate evaluation of the distribution functions f (x) 
and g(x) is required for todies not necessarily slender. For such 
todies a method is presented in references 3 and 9 for evaluating the 
upwash ty a stepwise procedure, and the stepwise procedure is extended 
to open-nose todies in reference 10. For a parabolic tody, a comparison 
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is presented in the EXAMPLES section of the incremental upwash due to 
angle of attack as computed by s lender-body theory (equation (15) with 
v =3 0 ) , by the more accurate stepwise procedure of reference 3, and by 
the intermediate procedure of using the s lender-body approximation for 
g(x) in the integral expression for the upwash. 

Lift due to interference upwash 

The calculation of the portion of the lift on a slender 
body due to the known upwash field w from another body is similar to 
the aft erbo dyjproblem in that the upwash varies along the body. The 
lift due to v is obtained in the same manner except that <f^ = 0. 

The upwash in each tra ns verse plane is assumed constant and equal to 
its value at the center line. The lift is then given by equation (10) 
with 

<J> ? = Ua R sin Q 
Ua ** 



according to s lender-body theory. The lift is thus 

L - 2qaAg — (27) 

Ua 

The same result is obtained for an open-nose body in a varying upwash 
field if the body is quasi- cylindrical (dR/dx - o) near the lip. 

The internal lift for an open-nose body of length substantially 
greater than p times the diameter with supersonic flow at the inlet 
is found from momentum consideration to be 


L = 2aqA 1 

provided there are no internal losses and the flow is discharged at 
the free-stream static pressure. Generally this will not be the case 
and terms must be included to correct for the internal axial forces. 

The lift acts essentially at the nose and the internal lift for an 
open-nose body in a varying upwash field is obtained simply by replacing 
the angle of attack by the ratio of the upwash to free-stream velocity 
at the nose. Thus the internal lift is 
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L = * 2a _ cji\.-T 
Ua 1 


(28) 


In this ana l ysis it is assumed that the portion of the lift of a 
body that results from the upwash generated by a second body due to 
the presence of the first body is negligible j in other words, only the 
primary interference effects are considered. In the section EXAMHES 
the secondary interference effects are exa min ed for a cylinder- cylinder 
combination. 

Buoyant lift due to interference 

The buoyant portion of the lift of a body due to the variable 
pressure field from another body is obtained from the second integral B 
of equation (5). Expanding the x-perturbation velocity u of the 
pressure field in a Taylor's series about the horizontal- plane through 
the center line of the body in question yields 



where Zq is the z-coordinate of the center line of the body. The 
terms in equation (29) involving even derivatives do not contribute to 
the lift, because they are symmetric about the z = z Q plane. If it 

is assumed that 


(z ~ zp ) 5 
31 



( z 



then that part of the u- velocity which contributes to the lift is 
approximately 


u *s (5 


< fb u\ 

- foL 

\ 'Z^Zq 


= R sin 0 




z=z r 


With the aid of the irrotationality condition ~ » -^2, the previous 


expression becomes 


bz bx 


u a R sin 0 



Z=sz, 


0 
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Inserting the value of u into the expression for the buoyant lift 
(term B in_equation (5)) and integrating ■with respect to 9 yield, 
■with w = w(x) , 



If the body is a cylinder or a quasi cylinder, 
to the pressure field is 

b = - a) 

\Ua Da, / 
or 

(c L ) -te(a.2L') 

\Ua Ua/ 


where A is the mean cross-sectional area. 

Equations (30 } and (3l) were obtained by expanding u in a 
Taylor 1 s series^ hence they are not valid in regions where a discon- 
tinuity in the u-velocity is present such as occurs, for example, at 
the intersection of the lip shock from an open-nose body with another 
body. The buoyant force, however, may be determined by applying equa- 
tion (30 ) or (3l), as the case may be, a suitable distance behind the 
intersection; then the buoyant force on the excluded area, associated 
with the shock, is evaluated separately and added thereto. The buoyant 
force associated with the shock is approximated herein by integrating 
the pressure as if it acted on the projected area (see following sketch) 
of the intersection of the shock with the body on the plane normal to 
the plane passing through the center lines of the bodies, and then 
resolving the force in the lift direction. 


the buoyant force due 

(31a) 

(31b) 
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It Is assumed that u and w are constant over the shaded surface 
area in the preceding sketch and have the values appropriate to the 
center line point P immediately behind the shock. This assumption 
implies that the body radius is relatively small compared with the 
separation of the bodies. For a slender body, the projected area is 
an ellipse with axes R and PR. By virtue of the assumption on u, 
the effective pressure acting on the area is constant, having the value 
given at the center line on the downstream side of the shock and zero 
on the upstream side. The value of the pressure on the downstream side 
of the Bhock is determined in appendix C using the same procedure as 
that used in reference 11 for the case of the zero angle of attack. 

The pressure coefficient is 


f - - f vt sia e ° (32) 

where the subscript 0 refers to the open-nose body. Since the projec- 
ted area on which this pressure acts is pjtR^, the buoyant force due to 
the discontinuity in u at the shock intersection with the body is 



sin^0Q 


(33) 
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Equation (33) together -with equation (30) or (3l) (-which apply beyond 
the shaded region in the sketch) yields the total buoyant lift when a 
discontinuity in the u- velocity is present. 

If the body is a cylinder or a quasi cylinder in the disturbed 
region of the pressure field so that equation (31b) may be used, the. 
expression for the buoyant lift with the lift due to the wave dis- 
continuity is given by 


2^ 2»j. AE A 

Wb - T'"u 73? 

where A is the projected area on which the shock pressure acts and 
w-j_ is the upwash immediately downstream of the wave. But for weak 

waves the flow deflection and pressure coefficient are related by 

Ap _ 2\ _ 2 ^ 

q p fi sin 0 q U 

where X is the flow deflection angle through the wave. Hence the 
contribution of the wave to the buoyant lift will be zero if 


= P sin e 0 (34) 

If the effective projected area is as assumed in the sketch, the rela- 
tion (34) is obtained. Thus, for this case, the correct buoyant lift 
is obtained by applying equation (3l) with w^ = 0. 

EXAMPLES 

Lift of Cone-Cylinder Body 

The lift of the cone part of the body as given in reference 3 is 



where x is the length of the cone. The lift of the cylindrical 
aftersection may be obtained from equation (19) with <p^ equal to the 

potential at the base_of the cone and cpg, the slender-body value 
(equation (14) with w = 0) . The potential on the body is 
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Substituting equation (14) with w *= 0 and equation (36) into equa- 
tion (19 ) and combining the result -with, equation (35) yield the lift 
of the cone- cylinder which, expressed in coefficient form based on the 
area of the cylinder, is 



In figure 3 the lift of the cone- cylinder and the lift of a cone 
are plotted as a function of the parameter x/f3R. As would be expected, 
the lift of the cone and the cone- cylinder approach the slender-body 
value as the cone angle decreases. 


Interference Lift Between Two Cylinders in Combination 

The lift of two semi-infinite cylinders of the same radius mounted 
in the vertical plane and starting at the same streamwise coordinate 
is computed to obtain some information concerning the magnitude and 
the method of computation of the interference lift between bodies. 

There is no loss in generality if the cylinders are assumed to be in 
the vertical, plane, and the computation is somewhat shorter. Also, 
if a symmetrically mounted vertical strut joins the bodies, the strut 
carries no lift and produces no significant Interference. 

The lift of one cylinder in the presence of the upwash field from 
the other cylinder is computed from equation (27). In figure 4 the 
upwash at various distances from a cylinder in the vertical plane as 
computed by the method of reference 10 Is presented. Rear the lip, 
values are improved by the method given in appendix C. For a short 
cylinder- cylinder combination, the appropriate value of the upwash to 
be used at the base would be the value given by figure 4. However, the 
upwash approaches the slender-body value rather rapidly downstream 6f 
the lip shock. In this example, therefore, the slender-body upwash 

Wg s* - UccR^/zq^ is used. The lift due to the upwash field in coef- 
ficient form and based on the cross-sectional area of the cylinder is 



3Y 


HACA 03ff 2669 


17 


°L = 



(38) 


where zq is the distance between the center lines of the cylinders. 

The interference lift given by equation (38) (which may be regarded 
as the first iteration of the total interference lift) can be considered 
to result from a certain distribution of doublets along the axis of the 
body. This doublet distribution produces an upwash field as though 
(in slender-body approximation) the upper body were moving in the nega- 
tive z-direction with velocity aU z^/lfi. The upwash induced at the 
lower body due to this motion at a distance Zq below the center line 
of the upper body is a fraction Zq2/r 2 Q f this. Thus the lift of the 

lower cylinder in the presence of the upwash from the upper cylinder 
due to the presence of the lower cylinder is given by 


°L 


= 2a 


R 


(39) 


The complete solution for the lift due to the upwash for both cylinders 
is then 



R c 


= - 4a 


1 + 


R 

V 


• + 


(-l) n+1 



(40) 


The v ali dity of this equation deteriorates as Zq diminishes because 

of the neglect of the gradient in upwash aeross each body. In the 
limiting case when the cylinders are touching (r/z^ a* l/2), equa- 
tion (40) is certa inly invalid quantitatively, but the behavior of the 
iterations is probably qualitatively correct. For this limiting case 
equation (40) becomes 


°L 




_ 1 _ 

64 


\ 4 

. . . ) a 

/ 5 


(41) 
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The magnitude of the terms of higher order than R 2 /zq 2 for the contig- 
uous case is 5 percent of the external lift of the cylinder alone. An 
exact -value for the upwash interference lift at small separation z Q> -within 

the frame-work of slender-body theory, may he found hy considering the 
two-d im ensional flow around two cylinders . This flow is equivalent to 
a doublet external to a cylinder. By use of the method of images 
(reference 12), the potential in series form (because each reflected 
doublet violates the boundary condition) may be obtained and hence 
the lift. The result for the interference lift is 



(42) 


For the contiguous case the interference lift due to the upwash field is 



= - 4a(l - 7^/12) 


+ (-l) n+1 

(n + 1) 


(43) 


The magnitude of higher-order terms for this method is about 7 percent 
of the external lift of the cylinder alone. The difference between 
equations (43) and (4l) is due to the fact that in obtaining equa- 
tion (40) the boundary condition is only approximately satisfied. 

A procedure analogous to that used to obtain equations (40) and 
(42) may be applied to obtain expressions for the buoyant lift. The 
buoyant lift is equal to the interference lift due to upwash. For this 
example, the internal interference lift is zero because the noses of 
the cylinders are ahead of the disturbance fields. 


The total lift of the combination, including the internal lift 
for the isolated cylinders, with terms of higher order t han R 2 /z Q 2 
neglected, is 



( 44 ) 
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The interference lift may "be expressed in terms of an interference ratio 


_ ^b dCL i 

da ~ da da R^ 

dCl h * 0 * 

da da 


(45) 


where is the lift of the isolated bodies. 


It is interesting to note that the interference lift ratio for the 
case of the cylinders in the horizontal plane would he 


i 



(46) 


Interference Lift Between Two Parabolic Bodies in Combination 

In this example, the interference lift between two identical 
parallel parabolic bodies mounted in the vertical plane and at incidence 
angle a with respect to the free-stream is computed. 

The upwash distribution about a parabolic body due to angle of 
attack may be computed either by slender-body theory, by the stepwise 
procedure of reference 3, or by the intermediate procedure of using the 
slender-body approximation for the g(x) in the integral expression 
for the upwash (equation (26)). The equation of a parabolic body in 
dimensionless form is 


R* » 


_ 2 _ 

F* 


(l - x*)x* 


(47) 


where 
R* pR/Z 

F* fineness ratio/p = F/p = 
x* x/Z 

Z length of body pointed at both ends 
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The upwash in the vertical plane through the body center line obtained 
by the intermediate procedure for the parabolic body defined by equa- 
tion (47) is 


_w_ 

Ua 



where Zq* = Pzq/Z. This equation is applicable between the Mach cones 
from the pointed nose and end of the body; that is. 


The upwash as computed by slender-body theory is 


w 

Ud 



(49) 


The comparison of equations (48) and (49) with the upwash computed 
using reference 3 for F* « 7-5 and 12.5 and Zq* = 0.075 and 0.125, 
respectively, is presented in figure 5. It is seen that the upwash 
given by equation (48) varies in the same manner with respect to magni- 
tude and position as that obtained using the method of reference 3. 

Since the computation with equation (48) is faster, equation (48) is 
used herein for the interference lift calculations due to the upwash 
from the interfering body. The sign of the interfering upwash acting 
at the base of the primary body is important because it determines 
whether or not the incremental lift due to interference upwash is 
favorable. The upwash given by slender-body theory will always have 
the same sign for a given plane of symmetry, whereas the upwash given by 
both of the more exact methods changes sign (see fig. 5). For the 
calculation of the buoyant lift, the simpler equation (49) (rather 
than equation (48)) has been used in equation (30) because, for inte- 
grated values, the difference between the various curves is unimportant. 


In figure 6 a contour plot of the interference ratio is presented 
for the case where the distance between the center lines of the bodies 
varies from Zq* = 0.075 to z* = 0.150 and one body moves downstream 

from x* = 0 to x* =0.3 with respect to the other parabolic body. 
The length of each body is 0.81. The maximum lift curve slope occurs 
when the body axes are farthest apart and one body is at the most 
rearward station of the positions considered. 



NACA TN 2669 


21 


Interference Lift Between Open- Nose Body and a 
Parabolic Body in Combination 

A contour plot of the interference ratio at a Mach number of 2 is 
presented in figure 7(a) for an open-nose body mounted above or below 
and downstream of a parabolic body. The bodies are parallel and at 
incidence angle a with respect to the free stream. All pertinent 
di m ensions are shown in the figure. In figure 7(b) another open-nose 
body identical to the first is symmetrically added to the configuration 
and the interference ratio is presented for this arrangement. 

For these plots , the upwash distribution about the parabolic body 
was computed from equation (48) with F equal to 15. For the open- nose 
body the upwash was determined Using the method of reference 10. For 
the computation of the buoyant force, equation (3l) was used with the 
mean radius determined by the amount of volume of the body in the dis- 
turbed field at a particular position. The lift of the isolated com- 
ponents was obtained using slender-body theory (equation (18) with 
w a 0). 

Figure 7(a) indicates that for the positions considered, the 
maximum lift-curve slope occurs when the open-nose body is at the 
farmost aft and outboard station from the parabolic body. The same 
result is noted in figure 7(b) for the parabolic body - two open-nose 
body configuration. 

As a matter of general interest, the interference lift at zero 
angle of attack is presented in figure 8 for the one open-nose con- 
figuration. When the open-nose body is mounted above the parabolic 
body, the lift is negative; and when it is mounted below, the lift has 
the same magnitude but the opposite sign. In the calculation of the 
lift the upwash about the parabolic body was determined from reference 8, 
and the upwash about the open-nose body was determined from reference 10. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, December 17, 1951 
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AFEEKDIS A 
SYMBOLS 


The 

A 

h 

B 

CL 

F 

^m-’^m 

f(x) 

S(x) 

i 

L 

l 

M 

n 

Ap 

5p 

1 

R 

r 

S 


following symbols are used in this report: 

c 

cross-sectional area of body of revolution, itR 
Tnayimrim semispan 

buoyant lift (second integral in equation (5)) 
lift coefficient, lift/qAg 

fineness ratio of parabolic body of revolution, 
length / maximu m thickness 

functions of x and r (defined in appendix B) 

local strength of source distribution 

derivative of local strength of doublet distribution 

interference lift ratio 

lift due to angle of attack and upwash (first integral in 
equation (5)) 

length of body pointed at both ends 

Mach number 

integer 

deviation of local pressure from free -stream pressure 

principal part of Ap contributing to lift 

free-stream dynamic pressure, i pU^ 

radius of body . 

maximum radius of body 

radial coordiante, 

wing area 
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t integration variable 

U free-stream velocity 

u,w perturbation velocities in x- and z-directions, respectively 

u,v,v perturbation velocity components due to interfering body 

V resultant velocity 

V Vg , V x components of resultant velocity in r-, 6-, and x- directions, 
respectively 


x,y,z 

a 

P 

r 

e 

£ 

p 

§ 

<p 

x 

* 


Cartesian coordiantes 
angle of attack 


cotangent of Mach angle of free stream, \ju Cj - 1 
adiabatic exponent 
tangential coordiante, sin _1 (z/r) 

integration variable (x-coordinate of sources or doublets) 
density 

total perturbation potential for primary body 
portion of $ contributing a nonzero value of lift 
portion of # corresponding to zero angle of attack 
portion of § due to sidewash 


Subscripts : 
0 


1 

2 

b 


coordinate measured from axis of body generating disturbance 
field 

start of body 
base of body 
bodies alone 
bodies in combination 
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i due to interference 

v -wing 

Subscript coordinates indicate partial differentiation irLtb respect 
to subscript variable. 

Primes denote ordinary differentiation. 
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APPENDIX B 

FOHM OF 'VELOCITY POTENTIAL FOR BODY OF REVOLUTION AND PRINCIPAL 

PARE OF PRESSURE CONTRIBUTING TO LIFT 

The total perturbation potential for the primary body is assumed 
to satisfy the governing equation for small-disturbance, isentropic 
irrotational flow 


_ * - S> ss o 

p xx v yy ^zz u 


(Bl) 


The general solution of equation (Bl) as given in reference 13 is 

CO OP 

$(x,u,e) = £ cos m Q F m (x,r) + sin m 9 C^fx^r) 




where 




m= 0 


f ffl (x - pr cosh t)dt 




cosh-1 

pr 


- r* 


esc 


x-Br 




0 ^(x - £) 2 - B 2 r 2 




m 0 


\ (B2) 


e-Fj* 


g^x - B^ cosh t)dt 


cosh _J_ 

Br 


= - i» 


HTf 


x-Br 




0 y(x - ?} 2 - B 2 r 2 J 


The linearized boundary condition on the body is 
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This boundary condition is correct to the order a cffi/cbc or (dR/dx) 2 . 
Because of the form of the boundary condition, the total potential 
must be given by 


§ = sin 0 + cos 0 F^_ + Fq 


“ «P + t + X 


(B4) 


(with v » v(x), w = w(x)) . The part of the pressure contributing to 
the lift will be obtained to the order of accuracy of the linear partial 
differential equation and the linearized boundary condition. The 
exact form of the pressure coefficient is given by the compressible 
Bernoulli equation 


AP 

CL 


2 

YM 2 


{[>• 


l^M 2 



JL } 

/ Y 2 ^! 

r-i 1 

\ U 2 /J 



(85) 


Expanding equation (B5) yields the following approximate expression 
for the pressure coefficient: 


AP 

1 



(B6) 


This equation is correct to the second order in 1 - V 2 /u 2 . The net 
velocity at any point is 


V 2 = Y x 2 + Y R 2 + Y@ 2 


where 


- 

y r » 


Yq =* 


U + u + § x 


(Ua + w)sin 0 + v 


(Ua + w) cos 0 + v 



(B7a) 

cos 0 + $ r 

(B7b) 

Q Q 

sin 0 + — 

(B7c) 


Substituting equation (B3) into (B7b) yields the net velocity in the 
radial direction of the body 


Y R = U 


dR 

dx 


(B8) 
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Substituting equations (B8), (B7a), and (B7c) into the first term in 
equation (B6) yields, on the body. 



The following terms in the approximate pressure coefficient yield 
nonzero values for lift when equation (B4) is substituted into (B9), 
the indicated differentiation is performed, the result is substituted 
into the lift integral 


Lift « 




— sin 0 R de dx 

q. 


(Bio) 


and the integration with respect to 


9 is performed: 


5£ = _ 2{ PxXx 2u 

q U U 2 U 2 u U 2 


(Bll) 


All other terms integrate to zero. If the next term in the 
of the pressure coefficient (equation (B6)) is included, by 
procedure 

8p _ 2 *s | [ 2 u | p 2 a 2 

q U U 2 TJ 2 U U 2 


expansion 
the same 


(bis) 


where the terms of higher order than those appearing in equation (Bll) 
have been neglected. It can be shown from slender-body theory that 
<Px/U, %/V) u/o, and X z /d 31,0 of order a dR/dx or (dR/dx) 2 . 


Neglecting terms of higher Older than a dR/dx or (dR/dx) 2 in 
equation (B12), where a and dR/dx are considered to be small, and 
of the same older, yields 


5p 2< ^x 2u 2u 2u ■ 

q 13 T" U U U 


(B13) 



28 


NACA TS 2669 


The component potentials >} f and X thus contribute only higher- 
order terms to the lift. _JEhe surviving camponent_in equation (B4) is <f } 
"Which is due to a and v alone; the sidewash v does not enter. 

Comparison of equations (B12) and (B13) indicates that at high 
Mach numbers some of the neglected terms may become as large as those 
retained. 


v 
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APPENDIX C 


EERTUERATION VELOCITY COMPONENTS ITT VICINITY OF LIP 


SHOCK FROM AN INCLINED OPEN-NOSE BODY 


The method used herein to estimate the velocity components in the 
vicinity of the lip shock from an open-nose "body at angle of attack is 
similar to the technique used in reference 11 for the linearized 
treatment of the lip shock from an open-nose "body at zero angle of attack. 
The boundary condition on the open-nose body in integral form may be 
written 



sin 0 


R 2 


x-pR 

f g'(S)(* - ?) 2 aj_ , 

0 Vo * - «) 2 - P 2 !* 2 


- Ua sin 0 


(Cl) 


where the lip is located at x = PR. In the vicinity of the lip, 

x - PR is very s mall compared with PR and thus equation (Cl) reduces 

to 


UaV2BR 

' p 2 


•f 


X-PR 


- PR - £ 


(C2) 


This is in the form of Abel's integral equation for the unknown function 
g'(x)j the solution, according to reference 14, Is 


g 'M («) 

when R is- assumed to be approximately constant in the vicinity of the 
lip. Introducing equation (C3) into the expression for the upwash 
(equation (26)), restricting consideration to the- vicinity of lip shock 
(x - pr very much less than pr), and integrating yield 

(C4) 

In the integration of the last term in equation (26), the function g(x) 
is given by the integral of equation (C3), which is 

g (x) = yj + o (c5) 

p^rt 

with the constant C equal to zero because g(o) 0. 
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The u-velocity component in the vici ni ty of the shock is obtained 
by inserting equation (C3) into the integral expression for the u- 
velocity and neglecting x - f3r compared -with (3r as before. The result is 

- ^ =■ - — sin 0 (C6) 

U q p Vr ' 1 

At the lip shock the upwash is zero when calculated by the method of 
reference 10. In figure 3 the upwash given by equation (C4) is joined 
by a dashed line to the values given by the method of reference 10 
downstream of the lip shock. The slope of the dashed line is about the 
same as that given by reference 10, thus indicating the validity of the 
method of this appendix. 
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Figure 1. - Body of revolution at angle of attaok in presence of disturbance velocity field 
I from another body. 
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Figure 3. - Comparison of lift of cone-cylinder witli lift 

of cone. 
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(a) pz Q /l, 0.075; F/p, 7.5. 


Method of reference 3 
Intermediate procedure, 
equation (48) 

Slender -body theory, 
equation (49) 




(b) Pz 0 /Z, 0.125; F/p, 12.5. 


Comparison of three methods of computing upwash in vei 
generated by a parabolic body at angle of attach. 
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Figure 8. - Contour plot of Interference lift coefficient at zero angle 
of attach: for open-nose body mounted above parabolic body. 
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